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ABSTRACT 



In a variability survey of M81 using the Large Binocular Telescope we have 
• discovered a peculiar eclipsing binary (My ~ —7.1) in the field of the dwarf 

^ ' galaxy Holmberg IX. It has a period of 271 days and the light curve is well-fit 

■ by an overcontact model in which both stars are overfiowing their Roche lobes. 
^ ■ It is composed of two yellow supergiants (V^ — / ~ 1 mag, Teg ~ 4800 K), rather 

than the far more common red or blue supergiants. Such systems must be rare. 
CN ■ While we failed to find any similar systems in the literature, we did, however 

\^ • note a second example. The SMC FO supergiant R47 is a bright {My ~ —7.5) 

cn 
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periodic variable whose All Sky Automated Survey (ASAS) light curve is well-fit 
as a contact binary with a 181 day period. We propose that these systems are 
the progenitors of supernovae like SN 2004et and SN 2006ov, which appeared to 
have yellow progenitors. The binary interactions (mass transfer, mass loss) limit 
the size of the supergiant to give it a higher surface temperature than an isolated 
star at the same core evolutionary stage. We also discuss the possibility of this 
variable being a long-period Cepheid. 

Subject headings: binaries: echpsing 



1. Introduction 



Although small in number, massive stars are critical to the formation and evolution 
of galaxies. They shape the ISM of galaxies through their strong winds and high UV 



fluxes, and are a raai or source of the heavy elements enriching the ISM (e.g, iMasseyl 12003 



Zinnecker fc Yor 



Ke 



20071. and referenc es therein). A large fraction of massive stars are found 



in binaries (e.g., iKiminki et al.l 120071 ). Eclipsing binaries are of particular use because they 
allow us to determine the masses and radii of the components and the distance to the system. 
Many young, massive eclipsing binaries have been found and studied in o ur Galaxy, the LMC 



and the SMC, prinaarily in OB associations and young star clus ters (e.g.. lBonanos et al.ll2004J : 
Peeples et al.ll2007l : iGonzalez et al.ll2005l : iHilditch et al.ll2005l ). The study of massive eclips- 
ing binaries beyond the Magellanic clouds has been limited until very recently, when vari- 
abihty searches using medium-sized telescopes with wide-field CCD cameras, coupled with 
spectroscopy using 8-met er class telescopes , have yiel ded the first systems with accurately 



measured masses in M31 ( Ribas et al. 2005 ) and M33 teonanos et al. 2006 ) 



We conducted a deep variability survey of M81 and its dwarf irregular companion, 
Holmberg IX, using the Large Binocular Camera (LBC) mounted on the Large Binocular 
Telescope (LBT), between January and October 2007. Holmberg IX is a young dwarf galaxy 
(age < 200 Myr), with a stellar populati on dominated by blue and red supergiants with no 
signs of old stars in the red giant branch (IMakarova et al.ll2002l ). Th e dwarf may have f ormed 
during a recent tidal interaction between M81 and NGC 2976 (e.g., Boyce et al. 200 J) . The 



gas-phase metal abundance of Holmberg IX of between 1/8 and 1/3 solar (e.g.. iMiller 



1995 



Makarova et al.l |2002| ) is consistent with this hypothesis (e.g., IWeilbacher et al.ll2003l ). A 
normal, isolat ed dwarf on the luminosity-metallicity relationship would have a metallicity of 
~ 1/20 solar Jhee et allbooel ). 



In this Letter, we report on the discovery of a 271 day period, evolved, massive eclipsing 
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binary in Holmberg IX using data from the LBT. The overcontact system is the brightest 
periodic variable discovered in our LBT variabihty survey. It has an out-of-eclipse magnitude 
of Vmax = 20.7 mag and is located at a = 09^57'^3r;14, 5 = +69°02'11" (J2000.0). In § 2 we 
discuss the observations and data reduction, in § 3 we present the light curve and the best-fit 
eclipsing binary model, and note a similar eclipsing binary in the SMC. In § 4 we discuss the 
results and their possible implications for Type II supernovae. Throughou t this Letter, we 
assume the HST Key Project distance to M81 of /i = 27.80 mag (3.6 Mpc; iFreedman et al. 



20011 ) as the distance to H olmberg IX, and cor rect only for a foreground Galactic extinction 
of E{B -V) = 0.08 mag JSchlegel et al.lll998h . 



2. Observations 

Holmberg IX was observed as part of a variability survey of the entire M 81 galaxy 



conducted between January and October 2007 with the LBT 8.4-meter telescope (|Hill et al. 
20061 ). using the LBC-Blue CCD camera jRagazzoni et allboodJciallongo et al.boOT^ during 



Science Demonstration Time. The survey cadence and depth (1 min single exposures, with 
> 3 consecutive exposures) are optimized to detect and follow-up Cepheid variables with 
periods between 10-100 days (V < 24 mag), getting better than 10% photometry in the B 
and V filters. We obtained 168 \^-band images on 24 different nights, and 87 5-band images 
on 13 nights. We coadded the 5-band images from each night (usually 3-4) to improve the 
signal-to-noise in the combined images. The seeing (FWHM) varied between 0'.'7 — 2'.'0 in V 
(median 1'.'4), and between I'.'O — 3'.'3 in 5-band (median 1'.'9). Our program did not request 
especially good image quality for these queue scheduled SDT observations. 

We also observed Holmberg IX as part of a variability survey of M81 conducted with 
the 8K Mosaic imager mounted on the MDM 2.4-meter telescope. The observations were 
obtained in 5 one-week runs between February 2006 and February 2007. All the images were 
obtained in \^-band using 15 min exposures. Due to weather loses and bad seeing, we ended 
up using only 36 images from 12 different nights. The typical seeing was ~ I'.'l. 



3. Light Curve 



We used the ISIS difference image analysis package (lAlardll2000l : iHartman et al.ll2004l ) 
to obtain the l^-band light curves of all the point sources detected in the LBT reference 
image. The detection of all point sources and the transformation of difference-flux light 
curves to instrumental magnitudes were done using the DAOPHOT/ALLSTAR package 
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( IStetsonl 119871 . Il992l ). After visual inspection of all the li ght curves of variable point-sources 
selected by standard criteria (rms and AoV significance; iHartman et al.l 120071 ). we detected 
~ 20 periodic variables in the field of Holmberg IX. These include Cepheids with periods 
of 10 — 60 days and one long-period variable. The analysis of the Cepheid PL-relation and 
the distance to Holmberg IX will be presented in a future paper (Prieto et al. 2008, in 
preparation). The brightest periodic variable is the pecuhar, long-period (P = 270.7 days) 
eclipsing binary we discuss here (hereafter VI). 

After discovering the binary in the LET data, we also ran ISIS and DAOPHOT/ALLSTAR 
on the MDM data to extract the light curve of the long-period binary. The variability data 
from LET and MDM were complemen ted with single-epoch archival im aging of the field ob- 
tained from the SDSS Data Release 6 (lAdelman-McCarthy et al.ll2007l ) in the gri bands (UT 
Nov. 30, 2003), and the HST/ACS Wide Field Camera (GO proposal 10605, PI E. Skill- 
man) in the F555W and F8UW filters (UT Mar. 23, 2006). The high resolution HST/ACS 
images (FWHM ~ O'.'l, corresponding to ~ 2 pc at the distance of Holmberg IX) show that 
the binary is spatially coincident with a stellar association in the dwarf galaxy. 

Figure [T] shows the phased V^-band light curve and B — V color curve of the eclipsing 
binary system. We include all the LET, MDM, SDSS and HST/ACS y-band photometry. 
The LET and MDM photometry have been calibrated using SDSS photometry of several 
relatively bright (r < 21.0 mag) and unsaturated stars in the field, transforr ning the qr 



magni tudes to standard BV magnitudes with the transformations presented in llvezic et al. 



(120071 ). The rms deviations of the absolute calibration are 0.02 — 0.03 mag for LBT-BV, and 
~ 0.05 mag for MDM-y. The SDSS gri photometry of the binary was extracted in the same 
way as for LET and MDM, using the DAOPHOT/ALLSTAR package to obtain instrumental 
magnitudes calibrated using absolute photometry of the bright stars in the field. Our g and 
r magnitudes of the binary from the SDSS data are 0.2 and 0.5 mag brighter, respectively, 
than the magnitudes reported in the SDSS-DR6 catalog, while the i-band magnitude agrees 
at the 1% level. We thin k this is due to prob lems in the SDSS photometry for faint sources 



in a crowded field (e.g., ISmolcic et al.l 120071 ). The details of the HST photometry can be 



found in Weisz et al. (2008, in preparation). 

In Figure [2] we show the position of the binary in the color magnitude diagrams (CMD), 
obtained from calibrated LET and HST/ACS photometry. The CMDs show the well- 
populated blue and red supergiant sequences in Holmberg IX. The binary is among the 
most luminous stars in this dwarf galaxy, with My ~ —7.1 mag, and it has clearly evolved 
from the main-sequence. With such a high intrinsic luminosity, the binary is bound to be 
massive. After correcting for Galactic foreground extinction, the B — V and V — / color s 
are consistent with an effective temperature of Tes = 4800 ± 150 K (IHoudashelt et al.ll2000l ). 
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Both components seem to be G-type yellow supergiants given the equal depths of the eclipses 
and the lack of color variations (see Figured]). 

We used the eclipsing binary model-fitting program NIGHTFALlH to model the V- 
band light curve. As shown in Figure [H we obtain a good fit to the light curve with an 
overcontact configuration where both stars are overflowing their Roche lobes. We fixed the 
effective temperature of the primary at Ti = 4800 K obtained from the colors. We assumed 
equal masses for the stars, a linear limb-darkening law, circular orbits, and synchronous 
rotation. We fit for four parameters: the Roche lobe filling factors, the inclination, and 
the temperature of the secondary. The time of the primary eclipse and the period were 
determined externally and were fixed for these fits. T he main paramet ers of the binary are 
listed in Table [H The light curve shows a hint of the lO'Connelll (119511) effect, in which the 
maxima (out-of-eclipse regions) show a difference in brightness (e.g., iPilecki et al.l 120071 ). 

We searched the literature for other examples of evolved, massive eclipsing binaries in the 
yellow supergiant phase and found non^. We also searched the available catalogs of eclipsing 
binaries in the LMC and SMC. The MACHO catalog of eclipsing binaries in the LMC 
(jPerekas et al.l 120071 ) contains 25 contact systems with red colors (i.e., evolved), (V ~ R) > 
0.5 mag, and periods > 200 days. However, these systems have absolute magnitudes My > 
—4 {V > 14.5 mag), that are ~ 3 mag fainter than the yellow supergiant eclipsing binary in 
Holmberg The All Sky Automated Survey (ASAS; Pojmanski 2002) contains complete 
Southern sky coverage for \^ < 15 mag. To our surprise, we found in the ASAS catalog a 
luminous (My ~ —7.5, Vmax ~ 11-5 mag), 181 day period contact eclipsing binary in the 
SMC. The star, SMC R47 {a = 01^29°^1?526, S = -72°43'20.2"), had be en spectrosc o picall y 
classified as an FO supergiant (Tcs ~ 7500 K) with emission lines by iHumphreyd (119831 ). 
The ASAS V^-band light curve of SMC R47, obtained between December 2002 and June 
2006, and the fit obtained with NIGHTFALL are shown in Figure O The best-fit eclipsing 
binary model requires a contact configuration, with non-zero eccentricity to account for the 
difference in timing between the eclipses. Even though we selected a clean part of the full 
ASAS light curve, there seems to be intrinsic variability from the binary components. The 
main parameters of this eclipsing binary are in Table [H While substantial ly hotter than the 
Holmberg IX binary, it does not lie on the SMC blue supergiant sequence ( iGrieve &: Madore 



^http : //www, hs .uiLi-hamburg.de/DE/Ins/Per/WlchmFiTin/Nightf all, html I 

^Note, however, that a possible Galactic counterpart is the contact binary BM Gas ( Fernie fc Evansl[l997 ). 
composed by an A7 lab supergiant (My ~ —6.3) and a late-type giant. 



Mennickent et al.l (|2006f) obtained spectroscopy of 17 "peculiar" periodic variables in the SMC from the 



OGLE database, and found a 184 day period eclipsing binary composed by two yellow supergiants (F5Ie 
G5-K0I). However, this system is 2 mag fainter than the binary in Holmberg IX. 
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1986). 



Discussion and Conclusions 



An eclipsing binary is the best explanation for the light curve of the brightest variable 
we have discovered in our LBT variability survey of the dwarf irregular companion of M81 
Holmberg IX. The other possible explanation for the periodic variability of VI is a long- 
period (P = 135 days) Cepheid. Such long-period Ceph eids (P > 100 days) have been 
observ ed in dwarf galax i es lik e the LMC and S MC (e.g.. iFreedrnan et al.l Il985h . NGC 55 



( e.g.. iPietrzvhs 



Pietrzvhski et al. 



ki et al.ll2006h . NGC 3 00 (e.g. 



Pietrzvhski et al. 



2002), NG C 6822 (e.g 



2004h . IC 1613 (e.g., kntonello et al.lll999[ ). and I Zw 18 jAloisi et al 



20071 ). The magnitude of VI is consitent with the magnitude of a Cephe id with P 
135 d ays {My — —7.0), extrapolating the period-luminosity relationship of iFouque et al. 
( 120071 ). However, while a few of these long period Cepheids have quasi-sinusoidal light 



curves that are nearly symmetric under a time reversal, they all have larger amplitudes 
in bluer bands {B amplitudes 1.3 — 1.6 times V amplitudes) due to the changes in the 



effect ive temperature as the star pulsates (e.g., iFreedman et al.lll985l : iMadore &: Freedman 
199ll ). Spectroscopy of VI, while challenging, would eliminate any remaining ambiguities in 
classifying this system. 

We can safely rule out the possibility that the eclipsing binary is in our Galaxy. Using 
the period, estimated effective temperature, magnitude at maximum, and assuming that the 
stars are of similar size in a contact configuration, we can use Kepler's law and their total 
surface brightness to e stimate a distance to the system of D ~ (Mtotaif^QY^^ Mpc, where 
Mtotai = M1 + M2 (e.g. . lGaposchkinlll962l ) . Conversely, we can estimate the total mass of the 
binary system by assuming the distance, Mtotai ^ 45 {D/3.Q Mpcf Mq. This should be taken 
as a rough estimate because of the overly simple model — to accurately constrain the total 
mass of the system, and its components, we need radial velocity measurements. Another 
piece of evidence that puts the binary system in Holmberg IX is its spatial coincidence with 
a stellar overdensity in the dwarf, observed in the HST/ACS images. 



We expected that such systems were rarqj, but were surprised to find none in the litera- 
ture. However, we found a similar eclipsing binary system in the SMC (SMC R47) searching 
through the ASAS catalog. From the absolute magnitudes of both binaries and their colors. 



*While the relative numbers of eclipsing binaries is a much more complicated problem, we note that 
the relative abundances o f red, blue and yellow supergiants is 4:13:1 for the Geneva evolutionary track 
(jLeieune fc Schaereiil2001r ) of a single, non-rotating star with M = 15 and Z = 0.004. 
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we estimate that at least one of the stars in each binary is ~ 15 — 20 MfD (main-sequence age 
~ 10 — 15 Myr) using the evolutionary tracks for single stars of iLejeune &: Schaererl (120011 ) 
(see Figure S]). 

The stellar evolutionary path of stars of a gi ven mass in bin ary systems can differ 
significantly from their evolution in isolation (e.g., iPaczyhskil Il97ll ) . In particular, binary 
interactions through mass loss, mass accretion, or common-envelope evolution , play a very 
important role in the pre-supernova evolution (e.g.. iPodsiadlowski et al.lll992l ). Most of the 
massive stars with masses 30 Mq ^ M > 8 M0 are expected to explode as supernova when 
they are in the red supergia nt stage, with a small contribu tion from blue supergiants (e.g., 
SN 1987A; IWest et al.l 119871 ). Surprisingly, iLi et al.l (120051 ) identified the progenitor of the 
Type IIP supernova 2004et in pre-explosion archival images and determined that it was a 
yellow supergiant with a main-sequence mass of ~ 15 Mq. Also, the posi tion in the CM P 
of the likely progenitor of the Type IIP supernova 2006ov (see Fig. 10 in iLi et al.l 120071 ) is 
remarkably similar to the position of the eclipsing binary in Holmberg IX (see Figure Hj). 

We propose that the binary we discovered in Holmberg IX and the binary found in the 
SMC@ are the kind of progenitor objects of supernovae like SN 2004et and SN 2006ov that 
appeared to be the explosions of yellow supergiants. A close binary provides a natural means 
of slowing the transition from blue to red, allowing the star to evolve and then explode as a 
yellow supergiant. As the more massive star evolves and expands, the Roche lobe limits the 
size of the star forcing it to have a surface temperature set by the uncoupled core luminosity 
and the size of the Roche lobe. It can expand further and have a cooler envelope only by 
becoming a common envelope system, which should only occur as the secondary evolves 
to fill its Roche lobe. This delayed temperature evolution allows the core to reach SN II 
conditions without a red envelope. 
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The authors thank the LET Science Demonstration Time (SDT) team for assembling and 
executing the SDT program. We also thank the LBC team and the LBTO staff for their 
kind assistance. We are grateful to Alceste Bonanos, Marc Pinsonneault, Lev Yungelson 
and Stan Woosley for helpful discussions and comments on earlier versions of this draft. We 
thank G. Pojmanski for making the excellent ASAS data publicly available. JLP and KZS 
acknowledge support from NSF through grant AST-0707982. 



possible earlier stage of thes e yellow supergiant binaries might be represented by the B-supergiant 
pair HD 1383 |Bovaiian et al-lbood ). 
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Fig. 1. — Phased V^-band light curve of the 270.7 day period echpsing binary in Holmberg IX 
{top panel) and its B — V color evolution {lower panel). The different symbols correspond to 
photometry from different telescopes: LBT {filled circles), MDM {open triangles), HST/ACS 
{filled square) and SDSS {open square). The solid line shows an overcontact eclipsing binary 
model that best fits the light V light curve. 
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Fig. 2. — Color Magnitude Diagrams of stars in the field of Holmberg IX obtained from the 
LBT-BV reference images {left panel) and HST/ACS-VJ single-epoch observations {right 
panel). The ACS CMD shows well-defined stellar sequences for the main sequence (MS) 
and the evolved blue (BSC) and red supergiants (RSG). The eclipsing binary {pentagon) lies 
between the blue and red supergiant sequences. 
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Fig. 3. — Phased V^-band light c urve of the long-p eriod, evolved eclipsing binary SMC R47 
obtained from the ASA S catalog (Poimanskill2002r) . The star was classified spectroscopically 



as an FO supergiant by iHumphreyd (1l983l ). The solid line shows the contact eclipsing binary 
model that best fits the light curve. The lower panel shows the residuals of the fit (observed- 
calculated). 
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Fig. 4. — CMD of Holmberg IX from HST/ACS V and / photometry. The connected filled 
symbols show the position of the evolved echpsing binaries in Holmberg IX {filled pentagons) 
and the SMC [filled triangles) at maximum and mini mum. The lines show ev olutionary 
tracks with extended mass-loss from the Geneva group (iLejeune fc Schaererll200ll ) for single 
stars with masses between 12-25 Mq, assuming two different metallicities: 1/3 solar (dashed) 
and 1/5 solar (solid). We use a distance modulus of = 27.80 mag to Holmberg IX, and 
Galactic color-excess E{B — V) = 0.08 mag, to put the evolutionary tracks in the diagram. 



- 15 - 



Table 1. Best-fit Binary Model Parameters. 



Parameter 


Holmberg IX VI 


SMC R47 


Period, P 


270.7 ±2.3 days 


181.58 ± 0.16 days 


Time of primary eclipse, Tp^jm 


2454186.0 ± 0.6 


2452073.1 ± 0.2 


Inclination, i 


55.7° ±0.6° 


82.2° ±0.2° 


Primary temperature, Ti 


4800 ± 150 K 


7500 ± 100 K 


Temperature ratio, T2/T1 


1.05 ±0.05 


1.17 ±0.02 


Eccentricity, e 


0.00 


0.039 ±0.002^ 


Roche Lobe Filling factors^ 


1.23 ±0.02 


1.02 ±0.02 


Semi-major axis, a'^ 


547 Rq 


418 Rq 



Note. — The mass ratio was fixed at g = 1 for fitting both fight curves. 

non-zero eccentricity is required to fit the difference in timing be- 
tween the primary and secondary ecfipses. The best-fit longitude of the 
periastron \s w — 168.6° ± 1.1°. 

^Ratio of stellar to Roche lobe polar radius for each star. 

■^Separation between the stars assuming a total mass of 30 M© for each 
system. 



